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An extension of the New Standard Model, by introducing a mixing of the low mass "active" neu- 
trinos with heavy ones, or by any model with lepton flavor violation, is considered. This leads to 
non-orthogonal neutrino production and detection states and to modifications of neutrino oscilla- 
tions in both, vacuum and matter. The possibility of the discovery of such effects in current and 
future neutrino oscillation experiments is discussed. First order approximation formulas for the 
fiavor transition probabilities in constant density matter, for all experimentally available channels, 
are given. Numerical calculations of flavor transition probabilities for two sets of New Physics pa- 
rameters describing a single "effective" heavy neutrino state, both satisfying present experimental 
constraints, have been performed. Two energy ranges and several baselines, assuming both the 
current (±2(t) and the expected in future (±3%) errors of the neutrino oscillation parameters are 
considered, keeping their present central values. It appears that the biggest potential of the discov- 
ery of the possible presence of any New Physics is pronounced in oscillation channels in which Ve, 
Ve are not involved at all, especially for two baselines, L = 3000 km and L = 7500 km, which for 
other reasons are also called "magic" for future Neutrino Factory experiments. 

PACS numbers: 13.15.+g, 14.60.Pq, 14.60. St 

Keywords: New Standard Model; Neutrino Oscillations; Neutrino Factory 



I. INTRODUCTION 



For several years, neutrinos are considered to be mas- 
sive particles and therefore the orthodox Standard 
Model (SM) with massless neutrinos must be extended. 
There exist two possibilities. Firstly, the extension of 
the SM can appear only at a very high energy scale, the 
GUT/Planck scale, and the non-zero neutrino masses are 
just the visible indication of such "high energy" physics 
in our "low energy" world. Such scenario is usually called 
the New Standard Model {vSM). Secondly, some New 
Physics (NP) can already be present at the TeV scale, 
that means at energies close to our present-day experi- 
mental facilities. The second of these possibilities is more 
appealing from both, the experimental, and the theoret- 
ical points of view. Such a "low energy" NP can par- 
ticipate in neutrino fiavor transitions, and so it could 
possibly be measured in future neutrino oscillation ex- 
periments. Then everything is in the hands of the pre- 
cision of the planned experiments. The bounds on the 
NP parameters, which arise from todays experiments, 
are too restrictive to give any good chance to see any 
effects in the present neutrino flavor transition data, tak- 
ing into account the fact that the present precision in 
determination of the neutrino oscillation parameters (of 
about 10% p]) effectively screens off any possible pres- 
ence of NP. However, the combined expected results from 
future neutrino facilities, like Beta Beam, Super Beam, 
and Neutrino Factory , should bring the neutrino oscilla- 
tion parameter errors down to about 1-3% [3| and there- 
fore give a chance for a discovery of effects, which possibly 
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could not be explained by the "present physics". 

The potential of the NP discovery is considered in this 
paper. There are many ways in which NP can mod- 
ify neutrino oscillations. The non-standard effects can 
directly change oscillation probabilities (the so called 
"damping signatures" or they can modify the oscil- 
lation amplitudes (by non-standard Hamiltonian effects 
[sj-flil), where both, oscillations in vacuum, and in mat- 
ter can be affected. These possibilities have extensively 
been examined in the existing literature. Thus, there are 
models with sterile 0| or/and heavy [H| neutrinos, gen- 
eral models with lepton flavor violation (LFV) 0|, non- 
standard interactions Q, flavor changing neutral currents 
general fermion interactions and mass varying 
neutrinos Next, there are models with non-unitary 
leptonic mixing violation of the Lorentz symmetry 
[131, violation of the principle of General Relativity 
and violation of the CPT symmetry • Finally, there 
are models which modify neutrino oscillations, i.e. models 
with neutrino wave packet decoherence neutrinos' 
decays (T3| and neutrino quantum decoherence (l8| . 

In this paper, we discuss one class of NP only, which 
can be obtained by mixing of the low mass "active" neu- 
trinos with heavy ones or by any model with LFV 
in both of which neutrinos interact with matter par- 
ticles by the left-handed (L-H) charge and neutral cur- 
rents only. In such models, the effective mixing matrix 
is non-unitary, resulting in non-orthogonal neutrino pro- 
duction and detection states. This non-orthogonality by 
itself modifles neutrino oscillations in vacuum. Apart 
from this, the neutrino interactions with matter particles 
are non-standard and so the oscillation effects in mat- 
ter are further modifled as well. Both these effects are 
here taken into account. Additionally, in our numeri- 
cal calculations, we assume both ±3% and ±2cr errors 
of today's vSM neutrino oscillation parameters, keeping 



2 



their present central values (see Chapter IIIip . The flavor 
transition probabilities were calculated for two energy 
ranges and several baselines [19]. Two of these baselines, 
L — 3000 km and L — 7500 km, are called "magic" for 
future Neutrino Factory experiments, as they are espe- 
cially useful for CP violation discovery {L — 3000 fcm) 
[20| or optimal for resolving the oscillation parameters 
degeneracy problem {L = 7500 km) . We have performed 
our numerical calculations of flavor oscillation probabili- 
ties for all available channels for two sets of NP param- 
eters which describe a single "effective" heavy neutrino 
state, both satisfying present experimental constraints. 
One of the easiest channels, from the experimental point 
of view is the i'f_i — > i^e one, but it will be difficult 
to observe any NP there (as it will also be in all another 
channels in which or are involved). In two other 
channels, — > and — > (and in the correspond- 
ing antineutrino channels) , the effects of NP are seen with 
the largest intensity, especially for the "magic" baselines 
mentioned above [22| . In the next Chapter we investigate 
how the NP modifies vSM oscillation transition prob- 
abilities. We discuss all leading terms, which give the 
new non-orthogonal production and detection states and 
which modify the neutrino coherent scattering on matter 
particles. We discuss also the reason why all channels 
with or Vr as the initial and final neutrinos (and the 
corresponding antineutrino channels), i.e. all channels in 
which I'e, t'e are not involved at all, are the most desired 
ones. Then, in Chapter ITTll we present the results of our 
numerical simulations. And finally, in Chapter ITVl we 
give our conclusions. In the Appendix [A] we collect all 
formulas for the flavor transition probabilities in constant 
density matter for all experimentally available channels. 



II. SEARCHING FOR NEW PHYSICS IN 
NEUTRINO OSCILLATION EXPERIMENTS 

If one takes into account that only relativistic neutri- 
nos are detected (and that only vector left-handed in- 
teractions are considered), the detection rate Nf^a of the 
neutrinos, coming from the produced Va neutrinos, 
factorizes into three parts, the production flux Na, the 



transition probability Pa-^p , 
tion Gf): 



and the detection cross sec- 



(1) 



Any physics beyond the i/SAI will modify all of these 
three parts (or even the above factorization will be made 
impossible). In this paper, we only discuss the modi- 
fications of the probability of neutrino fiavor transition 
from a production state | lyp) to a detection state I^d), 
leaving the modifications of the production and detection 
neutrino cross sections for future detailed considerations. 

In the frame of the vSM, the production or detection 
neutrino states are equal to the appropriate orthonor- 
mal neutrino flavor states \va)- The NP modifles this 



dependence and therefore the \vp,d) states are only ap- 
proximately equal to the l^a) states. Let us assume that 
the neutrinos are produced in the following process: 



e + X ^ ly + Y . 



(2) 



where £ is a charged lepton {£ — e,fi,T) and X,Y are 
hadrons. Then the normalized neutrino production state 
liyp) can be deflned as: 



j:t,A{e + x^u, + Y)\,,,) 



(3) 



where A{£ + X i^i + Y) is the amplitude for the process 
Eq.|[2|), in which the neutrino eigenmass state is pro- 
duced. The sum in Eq.|[3l) goes over all neutrinos with 
masses which are kinematically allowed. If particle 
spins were taken into account, instead of pure states of 
Eq.©, we would have to use mixed states described by 
an appropriate density matrix. 

Let us consider a NP model, in which besides three 
light i^SM neutrinos there are also heavier ones, which 
couple to the light charged leptons in a non-negligible way 
0. To be more precise, we assume the charge current 
Lagrangian in the following form {n > 3): 



^cc — 



E E 



2A/2sin6'vK " i 

7" (1 - 75) (U^)ai W- + h.c. 



(4) 



and similarly the neutral current Lagrangian in the fol- 
lowing form: 



C 



NC 



2sin(26l 



^ 7^ (1 - 75) VL^J Vj , (5) 



where 



(6) 



The n X n matrix XJ^, defines the mixing between the 
flavor and mass states. So, e.g. assuming three light and 
three heavy neutrinos we have: 



( ^1 \ 



(7) 



where all submatrices iJAy y and U ) have dimensions 
3x3. The three additional neutrino states, Ne,N^ and 
Nr, do not couple to charged leptons. For Majorana 
neutrinos the submatrices U and V, which explicitly en- 
ter the interaction Lagrangian Eq.|[4l), depend on 3n — 6 
moduli and 3n — 6 CP violating phases. For Dirac neu- 
trinos n — 1 phases can be eliminated, giving altogether 
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2n — 5 phases which, in principle, can enter the NP neu- 
trino flavor transition probabilities. If we assume that in 
the process given by Eq. ^ the energy conservation does 
not allow to produce heavy neutrinos Ni then, according 
to Eq.(l3]), the neutrino production state {\i^p)) is given 
by: 



1 



(8) 



Such states are normalized but not orthogonal. As the 
mixing of heavy neutrinos is small (IVa/p ^1), the ma- 
trix U is almost unitary. If we assume that a matrix U 
describes a unitary transition and is parameterized by 
the standard 3 mixing angles 6112 , 6113 , 6*23 and 1 standard 
Dirac CP breaking phase , then the orthonormal neu- 
trino flavor state Iva) is the following combination of neu- 
trino mass states It'i) 



(9) 



The U matrix can be parameterized by a matrix A close 
to the unit matrix 1: 



U^AU with A 1 - (5A , 



(10) 



and therefore the production state \iyp) in Eq.|[3l) is close 
to the eigenflavor state \i^a) and can also be decomposed 
in the orthonormal flavor basis: 



where the dia parameters are equal to: 

Ala 



Ell 



(11) 



(12) 



In general case, the values of parameters (1 — h)ia = 
{5K)pn = SXia depend on the production (detection) pro- 
cess P, [23I , and are bounded by the existing charged lep- 
ton data. The same parametrization as in Eq. lfTT]) was 
considered in Q , where the general lepton flavor violation 
NP model is discussed. In each row of the dia matrix, 
practically only one element has a non-negligible value, 
namely: 



\dee\ < 1 , and \dia \ « for a 7^ £ . 



(13) 



In general however, the 3x3 matrix SA can have all non- 
vanishing elements. Therefore, 9 moduH and 9 phases 
can generally parameterize any kind of NP. Not all phases 
play role in the transition probabilities. Five Majorana 
type phases do not enter any transition probability for- 
mula, hence only 4 phases remain. 

In the model considered here, the elements of the 6 A 
matrix are connected with the heavy neutrino mixing ma- 
trix V. From the unitary condition for the full matrix 
we get the following relation between A and V matrices: 



AAt 



(14) 



so, neglecting the SA6A^ term, we have: 

SA + SA^ ^ VV'' , 



or explicitly: 



6X 



S^*ja = Cap , where c^p = {VV'^) 



ap ■ 



(15) 



(16) 



Note that, these are not explicit formulas for individual 
5\ap elements. However, terms in form of the left hand 
sides of Eqs. lfT5|) . l|16p may often entirely describe the NP 
effects for modifled matter oscillations (see Eqs. (fT8l) . (|39|) 
below). If we need the knowledge of individual 5\ap el- 
ements, for example in order to calculate the new pro- 
duction and detection neutrino state modifications, then 
they need to be calculated from Eq. lfTOj) (see Eq. l(22|l be- 
low and the Appendix [A|. 

Up to now we have only considered the modifications 
of the initial and final neutrino states in the probabil- 
ity formula. Such modifications will change the neu- 
trino propagation even in vacuum. Yet, if neutrinos pass 
matter additional effects arise. The coherent neutrino 
scattering on matter particles is modified by NP and 
because of this: (i) neutrinos acquire different effective 
masses and (ii) their coherent scattering amplitude is 
modified. These effects can be parameterized by a NP 
effective Hamiltonian H^-^ (see [24], t25i]). The matrix 
representation of H^^ operator depends on the basis of 
states. Generally, in the basis of states of produced and 
detected neutrinos, the H^-^ operator is not represented 
by a hermitian matrix. However, it is represented by a 
hermitian matrix in the eigenmass basis {\vi)) and in any 
basis which is unitary transformed, e.g. in the basis of 
eigenfiavor orthonormal neutrino states given by Eq.|[9]). 
Therefore, in the orthonormal fiavor basis {\i'a)), we can 
write: 



H 



NP 



2E„ 



e 



Ett 



(17) 

where Ae — 2^/2GFNeEi, is the usual neutrino effective 
amplitude, which depends on the electron matter den- 
sity TVe, and Eap and Xap are NP parameters, moduli 
and phases, which describe the effective NP neutrino in- 
teraction with matter particles e, p and n. In the general 
case these parameters depend in a complicated way on 
the NP and matter properties. For uncharged, unpolar- 
ized and isotropic matter, parameters e are connected 
with c parameters in a simple way (see Eq. lfTO]) below). 
For example, in the frame of the model which we consider 
(Eqs.(l4]),|[5|), see also @), the effective NP Hamiltonian 
given by Eq. l|17p can be determined and, neglecting sec- 
ond order terms in 5A, is equal to: 



H^^ = ^{~A,[E{l)6A + 5A^E{l)] 



(18) 



-AnlSA + SA'^] 



where A^ = V^GpNnEi, depends on the neutron mat- 
ter density Nn (in the Earth matter An/A^ ~ 1/2) and 



4 



E{l)a/3 = SaeS/Se- Now, by Comparing equations Eq. (fT7|) 
and Eq. (fT8|) . we can find the following connection be- 
tween their parameters {(3 > a, Xaa = 0, it may often 
be the case that i5Ae^ = SXer = 0, and then Cafs entirely 
describe the NP effects for modified matter oscillations, 
compare Eq. l(39|) below): 

A 

EafSe^^"'^ = {-f- - SaeSf3e) Cq/3 - Saei^ " <^/3e) SXefS ■ (19) 

For high energy neutrino beams with w 0{GeV), 
the following two small parameters, which describe i/SM 
neutrino oscillations, are important (for lower energies, 
E,y « 0{MeV), a third small factor, Ae/Smli, would 
enter into the game, too): 



±0.03 and sin^ (26*13) < 0.05. 



(20) 



Also all Cq/3 and ea/3 parameters, which describe NP, 
are small. Therefore, we can expand neutrino oscillation 
probabilities in these small quantities, keeping only the 
leading first order terms. In this approximation, the full 
transition probability, for any fiavor and the baseline L, 
can be decomposed into two terms, the vSAI probability, 
and the correction to it given by NP: 



whereas the other terms satisfy: 

aUl3 



(27) 



E 

allp 



s prnass 



pint 

all 13 



In the next Chapter neutrino propagation in the Earth 
matter will be discussed. In our numerical calcula- 
tions of neutrino fiavor transitions we use the realistic 
PREM I [2HI Earth density profile model. However, ex- 
plicit analytical formulas for the fiavor transition proba- 
bilities can only be given for the case of constant density 
matter. Both NP corrections, which are important for 
neutrino transitions in matter, the 5PJj^^{L) and the 
|5-P™*^(i), are small, therefore their Hnear decomposi- 
tion in terms of a and sm(26'i3) is a very good approx- 
imation (in all formulas, we assume 5m\i = Sm^^i, and 
Srn^i = zLSrn^f^ + 6m1^i/2, where the upper/lower sign 
refers to the normal/inverted mass hierarchy (27l|): 



: pint _ pO 



sin(26'i3) B 



(28) 



and similarly 



Pp(a)^D{[j){L) ~ P^'^p{L) 



<;pNP 



[L). (21) 



X r)7nass 



sin(20i3)C 



/37 ■ 



(29) 



The NP correction probability we decompose again into 
two terms, the c term which is responsible for the initial 
and final neutrino state modifications (see Appendix [Aj 
and the e term which takes into account the NP infiuence 
on the coherent neutrino scattering in matter: 



5P, 



(22) 



The largest terms which are not suppressed neither by 
a nor by sin(20i3), namely the terms B^^ and C^^, do 
not appear in any i/^ nor i/g related channels. Such terms 



are only present in I'f^ 



and Vr Vr os- 



cillation channels (and in the corresponding antineutrino 
oscillation channels). For all such channels they are, up 
to the sign, the same and have the following form: 



The £ term consists of two terms also. The first 
one, which is responsible for an effective neutrino mass 
change, and the second one, which describes the addi- 
tional NP impact on coherent neutrino scattering with 
matter particles: 



(23) 



Generally, the production and detection states are not 
orthogonal: 



(24) 



and as a consequence the probability of neutrino oscilla- 
tion is not conserved: 



E Pp(a)~.P(!3) ^ 1 ■ 

alip 



(25) 



However, the neutrino oscillation probability of the vSM 
is normalized to one: 



all/3 



(26) 



B" = A, sin(4023){sin(2023) (e^M " ^rr) (30) 
+ 2cos(26'23) cos(xpr)£Mr}sin^(^) : 



and 



C° = AeAsin2(2023){-cos(2023)(£MM-£rr) (31) 
-f 2 sin(26'23) cos(x^t) e^r} sin(2A) , 



where 



A, 



:r-^ and A - 



'31 



2E^ 



(32) 



For the Vt- and transitions, we obtain: 

Blr = -Bl^ = and C^, - -C^^ - C° . (33) 

Unfortunately, these channels, wherein we can expect the 
largest NP effects, are not easily achieved experimentally. 
To see the NP effects in channels where they are sup- 
pressed, in the next Chapter, we discuss one of the easiest 
experimental channels, the v,, one. All non-leading 

terms in Eq. (|28|) and Eq. l|29p . together with all terms for 
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antineutrino and time reversal channels, are given in the 
Appendix [Aj 

In the model which we discuss, the c parameters are 
constrained from the existing experimental data (see 

Cee < 0.0054 , C,,^ < 0.0096 , Crr < 0.016 , 
ICe^l = Ic^el < 0.0001, ICerl - ICrel < 0.009, (34) 
Kr\ = ICr^l < 0.012. 

There are no constraints on the phases. The s parame- 
ters for the neutrino propagation in matter (Eq. lfT7|) ) are 
determined from relations given by Eq. (fT9|) . 

These are also the confinements which we use in the 
next Chapter, where we present some results of our nu- 
merical calculations of the differences of the transition 
probabilities between CP-conjugate channels: 

AP„%(L) = Po.^p{L) - P^_-p{L). (35) 

III. NEW PHYSICS IN FUTURE NEUTRINO 
OSCILLATION EXPERIMENTS 

In order to check the effects of the NP described above, 
the probability differences AP^^ (Eq. (|35|) ) for two en- 
ergy ranges and several baselines, have been calculated. 
Both, the energy ranges and the baselines have been cho- 
sen in prospect of the existing, planned, and feasible ex- 
periments. In view of Beta Beam and Super Beam ex- 
periments, the first energy range is Ey = 0.1 bGeV 
and the considered baselines are L = 130 , 295 , 810 km 
([29|). In view of Neutrino Factory experiments, the 
second energy range \s = 1 50 GeV and the con- 
sidered baselines are L = 732 , 3000, 7500 /cm (Q). In 
our numerical calculations of neutrino fiavor transitions, 
we use the realistic PREM I [2^ Earth density profile 
model, which assesses the actual matter density p and 
the actual electron fraction Y^, along the neutrino flight 
path in the Earth's interior. Then: 

A^AeV^] = 7.63xlO-M^][^][7?r;], 

g/cm-^ 0.5 GeV 

(36) 

A^[eV^] = 7.63xlO-M^^][l-ye][^]. 

g/cm-^ GeV 

I 

Two sets of V parameters, both satisfying present ex- 
perimental constraints given by Eq. l(34l) are discussed be- 



Note here that, for L < 874 km neutrinos pass only the 
first shell of the Earth's crust, with a constant density 
p = 2.Qg/cm^ and = 0.494, thus AeleV"^] = 1.96 x 
lO-^[E^/GeV] and A^[eV^] = 1.0 x IQ-'^lEjGeV]. 

The vSM oscillation parameters, together with their 
±2(7 errors (95% C.L., correlations among parameters are 
currently considered small), are taken from the current 
global best fit values (27| : 



sin2(0i3) = 0.9tgj X 10"2, 

Sml^i = 7.92 (1± 0.09) X 10"^ [el/^], 

sm\e^2) - 0.314 (lt^;i8), (37) 

5ml,^ = 2.4 (lt«-Ji)x 10-3 [ey^ 

sin2(^^23) = 0.44(ltl]-ji). 



In order to implement the effects of heavy neutrinos, 
as discussed in the previous Chapter, the matrix V must 
be introduced. As the number of heavy non-decoupling 
neutrinos is unknown, we parametrize it in a simplified 
way, using a single "effective" heavy neutrino state. In 
this way, the number of independent quantities parame- 
terizing the matrix Un in Eq.|[8l) are 6 moduli and 3 CP 
phases: 3 standard mixing angles (6'i2, ^13, ^23) + 1 stan- 
dard Dirac phase {Si^) + 3 new small NP mixing angles 
(014, 6*24, ^34) + 2 new NP Dirac phases ((524,(534). Then 



/ sin(0i4) \ 

V = cos(6'i4) sin(6'24) e-'^^^ , (38) 

\ 008(6*14) cos(6'24) sin(6'34) e-'^^' I 



and 



(39) 

I 

low: 

/ 0.001 \ / 0.01 \ 

(A) : 0.1 e-**^" , (B) : 0.01 e-^-^^^ . (40) 
yO.le-**^" / \ 0.1 e-'^^^ J 



1 - cos(6'i4) 

sin(6'i4) sin(6i24) e-**2* 1 - cos (6*24) 

sin(6'i4) cos(024) sin(6'34)e-^*^'' sin(6l24) sin(6'34) e^^^'^^''-''^'') 1-003(6134) 
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All calculations have been performed assuming the di- 
rect mass scheme only. As nothing is known about values 
of CP phases, we allow them to vary freely. 

For both V sets, we notice that, the biggest poten- 
tial of the discovery of the possible presence of any NP is 
pronounced in oscillation channels in which , are not 
involved at all, that is in t'^ Vt and I'^i (includ- 
ing the corresponding antineutrino channels). The effects 
are especially visible for two basehnes, L — 3000 km and 
L — 7500 km, which, for other reasons, are also consid- 
ered "magic" for future Neutrino Factory experiments. 
Moreover, comparing numerical results for these two sets 
of V parameters, we can clearly see that, as in Eq. l(40l B) 
the magnitude of the middle row is 10 times smaller than 
the corresponding magnitude in Eq. l(40l A). the NP effects 
for the second V set Eq. (|40l B) are smaller by a similar 
factor (compare Eq. l|4ip below), too. 

In order to find how the uncertainty of the estima- 
tion of the vSM oscillation parameters can mimic any 
possible NP effects, we have also performed calculations 
allowing all vSM oscillation parameters to vary by ±3% 
and ±2fT (note here that, the ±2cr test is also useful in 
qualitative estimation of the effect of the uncertainties in 
the Earth density profile on the vSM results). We have 
found that in general, in order to give a chance for a dis- 
covery of NP effects, it is required that vSM neutrino 
oscillation parameter errors should be diminished to the 
values expected in the future (of about ±3%). However, 
using neutrinos with energies i?^ > 15 GeV, together 
with the "magic" baselines mentioned above, it should 
even be possible to see NP effects with today's vSM neu- 
trino oscillation parameter uncertainties (of about ±2cr) . 
But, this chance depends on the actual magnitudes of 
the NP parameters and the actual precision of the ex- 
periment (which may, of course, be required to be much 
better than the allowed uSM neutrino oscillation param- 
eter errors). 

In Figdl the probability differences AP^]^^, for the 
second energy range, and the corresponding baselines set, 
are shown. The NP effects at L = 732 km, when scaled by 
a factor L/E^, can be used as an approximate estimation 
of the expected results for the first energy range with its 
considered baselines (for which we do not show pictures 
in this paper). 

In the left and middle columns of FigO the probabil- 
ity differences AP^fJ^^ for the second energy range and 
the two "magic" Neutrino Factory baselines, are shown. 
For this oscillation channel, the NP effect at L = 732 km 
is significantly smaller than in the channel at 

the same distance. This is also the case for this oscilla- 
tion channel in the first energy range with its considered 
baselines. Hence for this reason, we do not show the 
corresponding pictures in this paper either. 

Note that, in general, the transition probabilities of 
both, neutrino and antineutrino, oscillations depend on 
matter properties in a different way. However, at higher 
energies the vSM dependence is very much similar, thus 
any possible uncertainties in the Earth density profile 



cancel in the vSM probability differences (see the upper 
row of graphs, in both figures, for i?,^ > 15 GeV). In this 
way, any signal that the probability differences in this en- 
ergy range are distinctly different from zero will indicate 
that some NP exists. The statistical significance of such 
signal depends on its actual magnitude. For example, 
for the first V set Eq.lgOlA), and L = 7500 km, there 
exists a maximum at i?^ ~ 30 GeV, in which AP^J^^ w 
-AP^P^ « ±0.06, while P^^^ « P^^t^ ~ 0.5, and thus 
/^pCP I p ~ ±0.12 (see the middle row of graphs in both 
figures). The neutrino energy dependent magnitude of 
this effect can well be reproduced just by the largest 
term in the <5Pg^** NP correction (responsible for the 
effective neutrino mass change, see Eq. l|29p l. which is not 
suppressed neither by a nor by sin(20i3). Taking into ac- 
count the fact that P-a^'pi^iji ^e) = Pa^pi—Sij , —Ae), 
one can write (see Eq. (|3T|) ): 

AP^Sr « -AP^% « 2C° (41) 
« 4Ae sin^(26i23)A sin(2A) cos(x^r)epr. 

It should be noted, however, that the above formula does 
not reproduce the AP^J^^ well, in case one of the NP fea- 
sible CP phases (x^r — ^si — S24) is equal to tt/2 (see the 
dotted and dashed curves in Fig[T]). The reason is that 
in the above estimation we completely neglect several 
explicit CP asymmetry bracing terms, for example pro- 
portional to sin(xa^), which are relevant in the AP^_^^ 
case, but which never appear in the case of AP^fJ^^ (see 
the Appendix [A1 note also that, these asymmetry terms 
become dominant at short basehnes, and that is why the 
NP effect at L = 732 km is significantly smaller in the 
I'll channel than in the v-r channel at the 

same distance). From the above formula we can easily 
learn that, the magnitude of this NP effect is linearly pro- 
portional to both, the actual matter density through the 
term A^, and to the e^,- NP parameter (thus, as the value 
of the e^r parameter that results from the Eq. (|40l B) is 
10 times smaller than the one from the Eq. l(40l A). the 
NP effect for the second V set Eq. l|40I B) is smaller by a 
similar factor, too). 

As already mentioned, one of the easiest channels, from 
the experimental point of view, is the one, but it 

will be difficult to observe any NP there (as will be in all 
channels in which v^, or i^g are involved) . In the right col- 
umn of FigUl the probability differences AP^^g for the 
second energy range and the longest considered basehne, 
are shown (where the biggest effects are expected) . It can 
be seen that the NP effect is rather miserable in this oscil- 
lation channel, regardless of the V set given by Eq. l|40p . 
This conclusion holds also for all another similar oscil- 
lation channels and both energy ranges with the corre- 
sponding baselines. Moreover, these oscillation channels 
are sensitive to the (not so very well known) value of 
the sin^(0i3). In the other channels, without v^. and v^, 
the dependence on the sin^(6'i3) is small, giving a better 
chance to see the NP effects. 
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Figure 1: (Color online) The probability differences AP^J^^ for three baselines L = 732 , 3000 , 7500 km (each column of graphs 
corresponds to a single L), for the energy range E,, = 1 -{r 50 GeV. Calculations assuming the i/SM only are presented in 
the upper row of graphs, whereas results with NP are shown in the middle and bottom rows, except for the hashed band, 
identical in all rows, which corresponds to the current global best fit parameters with any feasible S13 value and uSM only 
(no NP). The dark (light) gray band in the upper row corresponds to ±3% (±2(t) deviations of the uSM neutrino oscillation 
parameters with any feasible S13 value. The light gray band in the middle (bottom) row corresponds to the current global 
best fit parameters with NP V parameters set Eq. l|40l A') (Eq. lfiOl Bll. with any feasible values of (5i3 , ^24 , <534 • Curves present 
in all graphs correspond to the current global best fit parameters with all feasible CP phases equal to (solid curves) and with 
exactly one of them equal to n/2 (dotted and dashed curves). 
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Figure 2: (Color online) The probability differences AP^^ for two baselines L — 3000 , 7500 fcm (the left and middle columns 
of graphs), and APl^Se for L = 7500 fcm (the right column of graphs) for the energy range Ei, = 1 50GeV^. Calculations 
assuming the vSM only are presented in the upper row of graphs, whereas results with NP are shown in the middle and 
bottom rows, except for the hashed band, identical in all rows, which corresponds to the current global best fit parameters 
with any feasible S13 value and uSM only (no NP). The dark (light) gray band in the upper row corresponds to ±3% (±2(t) 
deviations of the fSM neutrino oscillation parameters with any feasible (S13 value. The light gray band in the middle (bottom) 
row corresponds to the current global best fit parameters with NP V parameters set Ea. l|40I Al (Ea. l|40I Bll. with any feasible 
values of 5i3, ^24, <534. Curves present in all graphs correspond to the current global best fit parameters with all feasible CP 
phases equal to (solid curves) and with exactly one of them equal to tt/2 (dotted and dashed curves). 
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IV. CONCLUSIONS 

In the paper, by introducing a mixing of the low mass 
"active" neutrinos with heavy ones, we have investigated 
a possible New Physics (NP) scenario which is already 
present at the TeV scale, that means, at energies close 
to our present-day experimental facilities. In the pre- 
sented model (as also in any model with lepton flavor 
violation) , the effective mixing matrix is non-unitary, re- 
sulting in non-orthogonal neutrino production and de- 
tection states. This leads to the modification of the neu- 
trino oscillations in vacuum. Additionally, non-standard 
neutrino interactions with matter particles influence os- 
cillation effects also. First order approximation formu- 
las for the flavor transition probabilities, in constant 
density matter, for all experimentally available chan- 
nels, have been given. The possibilities of experimen- 
tal veriflcation of such model predictions have been dis- 
cussed in prospect of the existing, planned, and feasible 
Beta Beam, Super Beam, and Neutrino Factory exper- 
iments. Numerical calculations of flavor transition prob- 
abilities for two sets (satisfying present experimental con- 
straints) of NP parameters which describe a single "effec- 
tive" heavy neutrino state have been performed. They 
took into account two energy ranges and several base- 
lines, assuming both the current (±2(t) and the expected 
in future (±3%) errors of today's vSM neutrino oscilla- 
tion parameters, keeping unchanged their present central 
values. The realistic PREM I Earth density proflle model 
has been appHed. One of the easiest channels, from the 
experimental point of view, is the Vf^ Ve one, but it 
will be difficult to observe any NP there (as it will be 
also in all channels in which Ve or are involved). It 
appears that the biggest potential of the discovery of any 
possible presence of the NP is pronounced in oscillation 



channels in which v,,, vs are not involved at all, that 
is in Vt and t'^ (and in the corresponding 

antineutrino channels). The effects are especially visible 
for the two so called "magic" N eutrino Factory baselines, 
L — 3000 km and L — 7500 km. We have also found that 
in general, in order to give a chance for a discovery of 
the NP effects, it is required that i/SM neutrino oscilla- 
tion parameter errors should be diminished to the values 
expected in the future (of about ±3%). However, us- 
ing neutrinos with energies > 15 GeV together with 
the "magic" basehnes mentioned above, it should even 
be possible with today's uSM neutrino oscillation pa- 
rameter uncertainties (of about ±2(t). But, this chance 
depends on the actual magnitudes of the NP parameters 
and the actual precision of the experiment (which may, 
of course, be required to be much better than the allowed 
vSM neutrino oscillation parameter errors). 

Finally, it should be stressed that the full quantitative 
treatment of the NP effects in future facilities should be 
based on reaHstic observables related, for example, to the 
expected numbers of events at the corresponding facili- 
ties. However, as stated in the beginning of the Chap- 
ter mi this would require not only the knowledge of the 
NP generated transition probability modifications, but 
also the NP generated modifications of the production 
and detection neutrino cross sections need to be known. 
Some preliminary studies (see [^) suggest that the ex- 
pected effects on the last two terms can be of the same 
order as these on the first term (shown in this paper). 
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Appendix A: FIRST ORDER APPROXIMATIONS FOR FLAVOR TRANSITION PROBABILITIES 

Herein, we collect all formulas for the flavor transition probabilities, in constant density matter, for all experimentally 
available neutrino and antineutrino channels. As already mentioned, in all formulas, we assume Sm2i = Sm^^i, and 
Sm'^i = ±Sm^f^ + Sm'^^i/2, where the upper (lower) sign refers to the normal (inverted) mass hierarchy (27l |. 

Firstly, we show formulas required in order to calculate the c term in Eq. (|22l) . which is responsible for the initial 
and flnal neutrino state modifications (see also [l^l): 

SP^^piL) - {c^^ + C0p)P^%{L)-2ReM (Al) 
- 4 ^ i?e [{Sfy^^f,] sin2(^) - 2 J] /m [(<5f ^)L'>] sin(A,;,) , 

where 

(Sfy^p = -Y.(^^-'>U,^Upk + 5\f,,U^,U,k)KkU0. (A2) 

7 

- X! ('^'^"7 ^7* ^/3» + '^•^^7 ^"/s ^7i) ^f^'' ' 
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and 



2E^ 



(A3) 



Here, Uai and Smfj are, respectively, the effective mixing matrix elements and the effective mass square differences, 
coming from the diagonalization of the i/SM Hamiltonian in matter (without any NP, see, for example, |3l|), and 
dXaf3 can be calculated by putting SA = 1 — U W (compare Eq. lfTO|) ). 

Secondly, as described in ChapterlUl not all available channels were discussed in this paper. We were only interested 
in v^i~^ f ^ — > Vt and ~* channels, including the corresponding antineutrino transitions. In this Appendix, 
however, we collect all required equations. 

In any case, the total transition probability is decomposed into several terms according to 
Eqs. ([21]), ([221), {231), dsn),®. In general, the mass term i5P^^^ (Eqs.(l23]),ll29])) does not vanish only for chan- 
nels in which I'e, Ve do not enter. This term is also, up to the sign, the same for all such channels. Let us 
put: 



C" = A sin(2A) sin(26'i2) sin^(26'23) [003(6123) cos(xe^l.) £e^i - cos(xeT) sin(6'23)ee 



(A4) 



l-A_ 

For the transition, we obtain: 



A sin(2A) Sin^(26'23) [C0s((5i3 + Xe^) Sin(6'23) £ep + COs(6'23) C0s((5l3 + Xer) Eer] 



and two non-vanishing terms have the following form: 

1 



C0S(6'23) Sin(26'i2)<^ Eer cos(6'23) sin(6'23) 



COs(XeT)[^e C0S(2A) 



-Ae cos(2(Ae - 1)A) - 2{Ae - 2) sin2(leA)] + sin(xer)^e[sin(2A) 
- sin(2(l - Ae)A) - sm{2AeA)] 



+e 



efi 



2cos(xe;.)sin(^,A)[2(Ae - 1) 0082(6*23) sin(AeA) + ^e(sin((^e - 2)A) 



+ sin(AeA)) sin2(6'23)] + sin(xeM)^e sin2(6'23)[sin(2A) - sin(2(l - Ae)A) 
- sin(2leA)] 



(A5) 



(A6) 



(A7) 



Bf. 



Sin(6'23) \ Ser sin(6l23) C0S(6'23) 



. , ^. . , COs((5i3+Xer)[l+Ae 

{Ae ~ l)Al 

-{Ae ~ 1) cos(2A) - cos(2(Ae - 1)A) - (1 - Ae) cos(2leA)] 
+ sin(<5i3 + Xer)(le - l)[sin(2A) - sin(2(l - Ae)A) - sin(2leA)] 

+eeM cos((5i3 + XeMAe ~ 1) cos^ (6*23) (cos(2A) - cos(2(le " 1)A) 
+2&\T?{AeA)) + AAeSm^{{Ae - l)A)sin2(6l23)] 

+ sin(<5i3 + Xet.){Ae - 1) cos2(6l23)[sin(2(l - Ae)A) - sin(2A) + sin(2leA)] 



(A8) 



For the transition, all terms are different from zero. According to the previous discussion, the leading term 
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= (see Eq. l|33p ') and the two non-leading terms are given by: 



sin(2 6'23) 



Eep. sin(2 6'i2) sin(6'23) 



COs(Xe,0 



cos(2(-l + Ae) A) 



2 (-1+le) A, 

- cos(2 Ae A) (l~Ae + cos(2 623)) +2{2~Ae~2Al) 0082(6*23) sin2(A) 
+ {AAl -2Ae) sin2(6i23) sin2(A) - A^ sm{xe,,) [sin(2 A) + sin(2 (1 - A^) A) 
+ sin(2 Ae A) 



+eer cos(6'23)sin(2 6'i2) 



COs(Xer) 



cos(2 (-1 + Ae) A) - cos(2 A^ A) 



(-l + le + cos(2 6i23)) +2Ae cos2(6l23) sin2(A) -41^ cos(2 6l23) sin2(A) 



^2{Ae~2) sin2(A) sin2(023) 



+ Ae sin(xer) ^ sin(2 A) + sin(2 (1 - A^) A) 



f-sin(2v4e A) 

-4 (^-1 + Ae^ Al cos2(6ii2) cos(2 6I23) (2 A cos(A) - sin(A)) sin(A) 
sin(2 6*23) (e^^ - ffrr) + 2 cos(2 6*23) cos(x^^) e^r 



(A9) 



sin(2 9: 



23) 



Set, COS(023) 



2 cos(5i3 + Xep) sin(A) 



(AlO) 



2 + 

- (1 - 4Ae + 2l^) cos(2 6l23)j sin(A) - [A^ - cos(2 6'23)j sin(A - 2le A) 
-1 + Ae) sin(,5i3 + Xep) sin(2 A) - sin(2 (1 - A^) A) - sin(2 A^ A) 

2 cos((5i3 + Xer)sin(A) \{-A^ + {l-AA^ + 2Al) cos(2 6I23) ) sin(A) 



+eer Sin(6'23) 



{Ae + cos(2 6l23)) sin(A - 2 A A) - (^-l + A^ sm{5 



5l3 + Xe 



sin(2 A) - sin(2 (1 - A^) A) - sin(2 A^ A) 



For the mass correction terms, the first one C'^^ is given by Eq. l(33|) . and the two other universal terms are the 
following: 



(All) 



Also for the transition, all terms contribute. The modulus of the first interaction correction term is given by 
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Eq. l|33p . whereas the two remaining terms are as follows: 



_ . _ <^ £eii COs(XeM) COs(6'23) sin(2 6*12) ( -1 + Ae 

l + AA A, I LV 



1(2 1, A) cos4(023) + A, sin2(023) (-1 + cos(2 (-1 + le) A) sin2(023) 
-2 (-1 + A,) sin2(A) sin2(023)) + cos2(023) (l - ^ 



Sin^(023) 



+ (-1 + Ae) cos(2 (-1 + Ae) A) sin^(6'23) + A^ cos(2 A) 
2 sin2(A) sin2(6'23) + 2Ae sin2(A) sin2(6'23) + 21^ sin2(A) sin^ 



'23) + 2Ae sin^ 
cos(xeT) cos(6'23) sin(2 6'i2) 
2 cos(2 A) cos^ 



t'23 



1 

--£ 



'23) Sin(023) + 2 cos(2 (-1 + le) A) COS(023) sin3((?23) 
+4 cos(6l23) sin2(A) 810^(6*23) - sin(26'23) + ^ sin2(A) sin(46'2 
f (^-1 + cos2(6'i2) (2 A cos(A) - sin(A)) sin(A) sin(4f 

Sin(2 6*23) - Err) + 2 cos(2 ^23) cosixiir) e^r] 



(A12) 



-Bf,„ = 



(-1 + cos(2 A) cos2(023) + cos(2 (-1 + A^) A) Sin2(6'23)) 



ee^, COs((5i3 + Xep) 



2 ^ 



^^(-l + 2l.-cos(2M 



(A13) 



hC0S(^23) sin2(A) ((-l + le)Ae COs2(023) - (1 + ("3 + le) ^e) sin2(023)) 
in(2 6'23) + Eer C0s(5i3 + Xer) sin(6'23) 

-2 cos^(6'23) sin2(le A) sin(6'23) - 2 cos(6'23) sin2(A - A^ A) sin3(6'23) 
+ sin2(A) ((2 - Ae) le cos(2 623) + sin2(023)) sin(2 ^23)] | • 

The mass term for this channel has the opposite sign, in comparison to the previously discussed transition v^i ~^ Vt'- 

SPJT-^'^' = • (A 14) 

Finally, the i/SM probabilities are as follows: 



cos2(023) sin2(^. A) sin2(2 0i2) 

42 V 



(A15) 



sin2(2 0i3) 



-1 + A 
a sin(2 6'i3) 



(sin2(A-leA) sin2(023)) 

cos(5i3) cos(A) + sin((5i3) sin(A)] 



sin(Ae A) sin(A - A^ A) sin(2 6*12) sin(2 6*23) L 
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pSM 



sin^(A) sin^(2 6'23) - aA 003^(612) sin(2A) sin^(2 6i23) 

„2 



a 

4I2 



+^ <{ sin(A) sin(A - 2 A) sin2(2 6*12) sin2(2 6*23; 



sin2(2 0i3) 



sin(A) 2 ( - 1 ) A cos(A) 



4 (^-1 + A ^ 
+ sin( A) + sin( A - 2 A) sin^ (2 ^23 



+ - 



a sin(2 6'i3) 



COs((5i3) 



cos(2 6'23) sin(A) 



((-1 + 2^2^ sin(A) +sin(A- 2Ae A)) sin(2 0i2) sin(2 6'23) 



-2 sin(5i3) 



sin(A) sin(AeA) sin(A - A) sin(2 6'i2) sin(2 6'23) 



(A16) 



= l-sin2(A) sin2(2 023)+aAcos2(0i2) sin(2A) sin2(2 023) 



sin" (2 013) 
4 (-1 + A, 



-4 sin"(A- A) sin'^(6'23) 



siir(A) + [-l + Ae] A sin(2 A) - sin^(yle A) sin'' (2 6*23) 



sin2(2 6li2) <^ 4 cos'^(6'23) sin"(Ae A) + sin^CA - Ae A) sin2(2 6'23; 



(A17) 



4A2 

^ a£in(2M_ sin(2 0i2) (sin2(A) (l" cos(2 ^^23) + sin2((?23)) sin(2 023) 

- sin(023) (2 cos3(023) sin2(l. A) + sin2(A - A^ A) sin(023) sin(2 ^23; 

The formulas for the time reversed channels, i^e and Vr ^ v^i^ are obtained from the above formulas after the 

replacements A ^ — A and A^^ —^ik- Finally, if no scalar nor pseudoscalar neutrino interactions are considered 
(like in the model which we present here, see also for antineutrinos the interaction Hamiltonian H is replaced 



by —H*, thus the transitions formulas for channels vji 



Vf and Vfi vp, are easily obtained after the 



replacements S13 — > — (5i3, A^ ^A^, and SXap ^^a0 (hence also Xap — > ~Xaf3)- This completes the set of i/SM 
transition probability formulas and the NP corrections to th em, for all experimentally available channels. 
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